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Abstract 

By starting from the stochastic Schrodinger equation and quantum 
trajectory theory, we introduce memory effects by considering stochastic 
adapted coefficients. As an example of a natural non-Markovian ex- 
tension of the theory of white noise quantum trajectories we use an 
Ornstein-Uhlenbeck coloured noise as the output driving process. Un- 
der certain conditions a random Hamiltonian evolution is recovered. 
Moreover, we show that our non-Markovian stochastic Schrodinger 
equations unravel some master equations with memory kernels. 
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1 Introduction 



The problem of how to describe the reduced dynamics of a quantum open 
system interacting with an environment is very important [TH3]. More and 
more apphcations demand to treat dissipative effects, tendency to equihb- 
rium, decoherence,... or how to have more equihbrium states, survival of 
coherences and entanglement,... in spite of the interaction with the ex- 
ternal environment. The open system dynamics is often described in terms 
of quantum master equations which give the time evolution of the density 
matrix of the small system. When the Markov approximation is good (no 
memory effects) the situation is well understood: if the generator of the 
dynamics has the "Lindblad structure" , then the dynamics sends statistical 
operators into statistical operators and it is completely positive |4|5] . 

However, for many new applications the Markovian approximation is not 
applicable. Such a situation appears in several concrete physical models: 
strong coupled systems, entanglement and correlation in the initial state, 
finite reservoirs... This gives rise to the theory of non-Markovian quantum 
dynamics, for which does not exist a general theory, but many different 
approaches [6H23]. 

Non Markovian reduced dynamics are usually obtained from the total 
dynamics of system plus bath by projection operator techniques such as 
Nakajima-Zwanzig operator technique, time-convolutionless operator tech- 
nique |2|10j . correlated projection operator or Lindblad rate equations [TH 
119]... These techniques yield in principle exact master equations for the 
evolution of the subsystem. For example Nakajima-Zwanzig technique gives 
rise to an integro-differential equation with a memory kernel involving a re- 
tarded time integration over the history of the small system. However, in 
most of the cases the exact evolutions remain of formal interest: no analytic 
expression of the solution, impossible to simulate... Usually, some approx- 
imations have to be used to obtain a manageable description. But as soon 
as an approximation is done, the resulting equation can violate the complete 
positivity property; let us stress that the complete positivity (and even pos- 
itivity) is a major question in non-Markovian systems |18|22j . 

The easiest way to preserve complete positivity is to introduce approx- 
imate or phenomenological equations at the Hilbert space level, an approach 
which is useful also for numerical simulations. We can say that in this way 
one is developing a non-Markovian theory of unravelling and of "Quantum 
Monte Carlo methods" |2|6H9|12|24] . In the Markovian case such an ap- 
proach is related to the so called stochastic Schrodinger equation, quantum 
trajectory theory, measurements in continuous time. It provides wide ap- 
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plications for optical quantum systems and description of experiments such 
as photo- detection or heterodyne/homodyne detection |24fl31] . In the non- 
Markovian case, an active line of research concentrates on a similar inter- 
pretation of non-Markovian unravelling. In this context, the question is 
more involved (for example, when complete positivity is violated the an- 
swer is hopeless) and remains an open problem. For the usual scheme of 
indirect quantum measurement it has been shown that in general such an 
interpretation is not accessible |21|23] . Actually, only few positive answers 
for very special cases have been found and this question is still highly de- 
bated |12ll3l2UI21j . 

Our aim is to introduce memory at Hilbert space level, in order to guar- 
antee at the end a completely positive dynamics, and to maintain the pos- 
sibility of the measurement interpretation. Our approach is based on the 
introduction of stochastic coefficients depending on the past history and on 
the use of coloured noises [32H34] . 

In Sect. [2] we introduce a special case of stochastic Schrodinger equation 
with memory. The starting point is a generalisation of the usual theory of 
the linear stochastic Schrodinger equation [3TH33] . based upon the intro- 
duction of random coefficients. This approach introduces memory effects 
in the underlying dynamics. The main interest is that the complete posit- 
ivity is preserved and a measurement interpretation can be developed. We 
present this theory only in the context of the diffusive stochastic Schrodinger 
equation. 

In Sect. Owe attach the problem of memory by introducing an example 
of coloured bath and we show that we obtain a model of random Hamiltonian 
evolution [34], while we remain in the general framework of Sect. [21 

Finally, in Sect. [H by using Nakajima-Zwanzig projection techniques, 
we show that the mean states satisfy closed master equations with memory 
kernels, which automatically preserve complete positivity. Moreover, we can 
say that the stochastic Schrodinger equations of the previous sections are 
unravellings of these memory master equations. 

2 A non Markovian stochastic Schrodinger equa- 
tion 

The linear stochastic Schrodinger equation (ISSE) is the starting point to 
construct unravelling of master equations and models of measurements in 
continuous time [27|31j . By introducing random coefficients in such equa- 
tion, but maintaining its structure, we get memory in the dynamical equa- 
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tions, while complete positivity of the dynamical maps and the continuous 
measurement interpretation are preserved |32|33j . To simplify the theory we 
consider only diffusive contributions and bounded operators. 

Assumption 1 (The linear stochastic Schrodinger equation). Let "K he a 
complex, separable Hilbert space, the space of the quantum system, and 
(il,3", (3'j),(Q) be a stochastic basis satisfying the usual hypotheses, where 
a d-dimensional continuous Wiener process is defined; Q will play the role 
of a reference probability measure. The ISSE we consider is 



dip{t) = K{t)iP{t)dt + Rj{t)iP{t)dWj{t), 



Let us denote by T(!K) the trace class on [K, by S(?{) the subset of the 
statistical operators and by L{!K) the space of the linear bounded operators. 

Assumption 2 (The random coefficients). The coefficient in the drift has 
the structure 



i^(t) = -i//(t)-i^i2,(t)*i?,(t). 



(2) 



The coefficients H{t), Rj{t) are random bounded operators with H{t) = 
H{t)*, say predictable caglad processes in (^^,3", (3"*), Q)- 
Moreover, VT > 0, we have 

rT 

EQ[||F(t)||]dt<+oo, (3a) 



{2E/ 11^. 



exp 



< +00. 



(3b) 



Theorem 1. Under AssumptionsUll^ the ISSE ([T]) has a pathwise unique 



solution. The square norm 
by 



iJo\\ exp 



E 



7' 





mj{s)dWj{s) - ^ I mj{s)'^ds 



is a continuous positive martingale given 
1 /•* 

(4) 

(5) 
(6) 



mj{t) := 2Re(^ip{t)\Rj{t)i;{t)^ 
\v (fixed unit vector), if ||V'(t)|| = 0. 
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Proof. Assumptions [H [2] imply the Hypotheses of [32, Proposition 2.1 and 
Theorem 2.4], but Hypothesis 2. 3. A of page 295. According to the discus- 
sion at the end of p. 297, this last hypothesis can be substituted by ()3bp . 
which implies Novikov condition, a sufficient condition for an exponential 
supermartingale to be a martingale. Then, all the statements hold. □ 

Remark 1. By expression (j4]) we get that on the set {HV'oll > 0} we have 
||V'(t)|| > Q-a.s. This means that, if ipQ ^ (Q-a.s., then the process ilj{t) 
([6]) is almost surely defined by the normalisation of ^l^{t) and the arbitrary 
vector V does not play any role with probability one. 

Remark 2. Let us define the positive, 'J(J{)-valued process 

a{t):=\m){m\- (7) 

By applying the ltd formula to {^p(t)\a^p{t)) , a G 'C{!K), we get the weak- 
sense linear stochastic master equation (ISME) 

d 

da{t) = C{t)[a{t)]dt + '^nj{t)[a{t)]dWj{t), (8) 
i=i 

nj{t)[p]:=R,{t)p + pR,itr, (9) 

Cit)[p] = -i[H{t),p]+Y.(^j(^)pRj(^y -URj{trR,{t),p}] ; (10) 

i=i ^ ^ 
C{t) is the random Liouville operator [32, Proposition 3.4]. 

Assumption 3 (The initial condition). Let us assume that the initial con- 
dition ijjQ is normalised, in the sense that Eq HV'oIP = 1- Then, go := 
Eq [|'i/'o)(V'o|] G S(^K) represents the initial statistical operator. 

Remark 3. Under the previous assumptions p{t) := ||'i/'(^)||^ is a positive, 
mean-one martingale and, VT > 0, we can define the new probability law 
on (fl, Jt) 

VF G Jt P^{F) := Eq\p{T)lF]. (11) 

By the martingale property these new probabilities are consistent in the 
sense that, for < s < t and F G J^, we have P*(F) = P^(F). 

The new probabilities are interpreted as the physical ones, the law of 
the output of the time continuous measurement. Let us stress that it is 
possible to express the physical probabilities in agreement with the axiomatic 
formulation of quantum mechanics by introducing positive operator valued 
measures and completely positive instruments [32ji33j . 
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Remark 4. By Girsanov theorem, the d-dimensional process 



Wj{t) ■.= Wj{t) 



mj{s)ds, j = l,...,d, t£[0,T], 



(12) 



is a standard Wiener process under the physical probabihty P"^ \32\ Propos- 
ition 2.5, Remark 2.6]. 

By adding further sufficient conditions two more important equations 
can be obtained. 



Assumption 4 ( |32^ Hypotheses 2.3.A]). Let us assume that we have 



sup / Rj{s,uj)* Rj{s,uj) 
weQ Jo " 



ds < +00. 



(13) 



Theorem 2. Let Assumptions Ui]^ hold. Under the physical probability 
the normalized state ip{t), introduced in Eq. ([6]), satisfies the non-linear 
stochastic Schrodinger equation (SSE) 

d^{t) = \Rj{t) - Renj (t,^(t)) i^{t)dWj{t) + K{t)^l){t)dt 



Y (Renj (t,Vi(t))) Rj{t) - ^ (Renj (t,V'(t) 



mdt, (14) 



where njit, x) := {x\Rj{t)x) , Vt G [0, +oo), j = 1, . . . ,d, x £ !K. 
Moreover, the process ( a priori or average states ) defined by 

r?(t) := EQ[(7(i)], or Ti {ar^{t)} = ^^[{^P{t)\ail:{t))] Va G £(J{) (15) 

satisfies the master equation 



Vit) = Qo 



Teq [C{s)[a{s)]]ds. 
Jo 



(16) 



Proof. One can check that all the hypotheses of [32| Theorem 2.7] hold. 
Then, the SSE for ip{t) follows. 

As in the proof of |32l Propositions 3.2], one can prove that the stochastic 
integral in the ISME ([8]) has zero mean value. Then, Eq. ()16p follows. □ 

Note that, by the definition of the physical probabilities, we have also 
Tr {ar/(t)} = EpT [(V'(t)|a?/'(t))l , Va G £(:K), Vt, T : < t < T. (17) 



6 



The SSE ()14p is the starting point for numerical simulations; the key point 
is that norm of its solution ijj{t) is constantly equal to one. We underline 
that Eq. (jl6p is not a closed equation for the mean state of the system. In 
the last section we shall see how to obtain, in principle, a closed equation 
for the a priori states of the quantum system. 



The finite dimensional case 



If we assume a finite dimensional Hilbert space and we strengthen the con- 
ditions on the coefficients, we obtain a more rich theory. We discuss here 
below the situation [33j . 

Assumption 5. 1. The Hilbert space of the quantum system is finite 
dimensional, say !K := C". We write M„(C) for the space of the linear 
operators on into itself (n x n complex matrices). 

2. The coefficient processes Rj and H are M„(C)-valued and progressive 
with respect to the reference filtration. 

3. The coefficients satisfy the following conditions: for every T > there 
exist two positive constants M{T) and L{T) such that 



sup sup 

wente[o,T] 



< L(r) < oo. 



(18) 



sup sup \\H{t,uj)\\ < M{T) < oo. 

wGH ie[0,T] 



Under Assumption [5] it is possible to prove existence and pathwise unique- 
ness of the solution of the ISSE just modifying classical results for existence 
and uniqueness of the solution for stochastic differential equation with de- 
terministic coefficients. 

Moreover, it is possible to prove that the solution of ISSE fulfil some 
LP estimate: in this point the finite dimension of the Hilbert space plays a 
fundamental role because the bounds we obtain for the process ip{t) involve 
constants depending on n. 

Obviously, in this context the martingale property of the norm of the 
solution is still valid and so one can define the consistent family of phys- 
ical probabilities. It is also possible to introduce the propagator of the 
ISSE, that is the two times M„(C)-valued stochastic process A{t,s) such 
that A(t,s)^p{s) = tp{t), for all t,s > s.t. s < t. We are able to ob- 
tain a stochastic differential equation (with pathwise unique solution) for 
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the propagator and, by means of it, to prove that the propagator takes 
almost surely values in the space of the invertible matrices. The es- 
timates for ip{'^) 'ire useful to obtain estimates on A{t,s). Further- 
more, the propagator satisfies the typical composition law of an evolution: 
A{t, s) = A{t, r)A{r, s) for ah t, r, s > s.t. s <r <t. 

The almost sure invertibility of the propagator guarantees that the pro- 
cess il){t) can be almost surely defined and that this process satisfies, under 
the physical probabilities, a non linear SSE, similar to Eq. ([HI). It is possible 
to prove that in this case the SSE has a pathwise unique solution. 

When we go on extending the theory to the space of the statistical op- 
erators, we can take as initial condition a random statistical operator or a 
deterministic one. We define the process a{t) as in Eq. ([7]) and, by using 
the Ito formula, we obtain an equation formally similar to the ISME, but 
in this case we are able to prove the uniqueness of its solution given the 
initial statistical operator qq (the existence comes out by construction). In 
this way we can say that the ISME is the evolution equation of the quantum 
system, when the initial condition is a deterministic (or even random) stat- 
istical operator qq. We can introduce the propagator of the ISME, which 
is a two times-linear map valued stochastic process, say A(t, s), such that 
A(t, 0)[£>o] = cr(t), A(t, s) = A{t, r) o A(r, s), for all t, r, s > s.t. s <r <t 
and A(t, s)[r] = A{t, s)TA{t, s)* . From the last expression of the propagator 
of the ISME, it comes out that this is a completely positive map valued 
process. 

Also in this case we can introduce a consistent family of physical probab- 
ilities. Indeed, the process Tr{a{t)} is an exponential mean-one martingale 
that can be used to define the new probability laws, as we did in the Hilbert 
space. 

It is then possible to define the normalisation of a{t) with respect to its 
trace, 

and, under the physical probabilities, we have the following non linear equa- 
tion for g{t), with pathwise unique solution 

d 

d0{t) = c{tMt)]dt + Y,{T^j{tMt)]-vj{t)e{t)}dw,{t), t>o 

i=i 

£•(0) = £'0 , 

(19) 
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where Vj{t) := Tr{{Rj (t) + Rj (t)*) g{t)} , and W{t) is a Wiener process under 
the physical probabihties defined by 




3 Random Hamiltonian 

In the previous section, we have presented a non Markovian generahsation of 
the usual difi'usive ISSE by using random coefficients to introduce memory. 
In this section we adopt an alternative strategy and we start with a usual 
ISSE with non random coefficients, but driven by a coloured noise; in this 
way the memory is encoded in the driving noise of the ISSE, not in the 
coefficients. As we shall see, this model too turns out to be a particular case 
of the general theory presented in Section [2j Moreover, the new ISSE will 
be norm-preserving and will represent a quantum system evolving under a 
random Hamiltonian dynamics, while the Hamiltonian is very singular and 
produces dissipation. 

As our aim is just to explore some possibility, we keep things simple and 
we consider a one-dimensional driving noise and two non-random, bounded 
operators A and B on'K in the drift and in the diffusive terms. The starting 
point is then the basic linear stochastic Schrodinger equation 



The simplest choice of a coloured noise is the stationary Ornstein-Uhlen- 
beck process defined by 



where {W{t)) is a one dimensional Wiener process, defined on the stochastic 
basis {Q, 3", 3"t, (Q) and Z is an 3"o-measurable, normal random variable with 
mean and variance 1/(27). The Ornstein-Uhlenbeck process {X(t)) is a 
Gaussian process with zero mean and correlation function 



dtp{t) = A^{t)dt + B^{t)dX{t). 



(20) 




7 > 0. 



E^[X{t)X{s)] 



27 



(21) 



it satisfies the stochastic differential equation 



dX{t) 



-/X{t)dt + dW{t) 



X{0) = Z. 



(22) 
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Formally, Eq. ()20p is driven by the derivative of the Ornstein-Uhlenbeck 
process (heuristically, dX{t) = X{t)dt), whose two-time correlation is no 
more a delta, as in the case of white noise, but it is formally given by 
Eq[X{t)X{s)] = 5{t - s) - ^e-'^l*-'^!. Note that the Markovian regime is 
recovered in the limit 7^0. 

It is then straightforward that Eq. ()20p can be rewritten in the form 

d^{t) = {A- 7X(t)S)V(t)dt + BiP{t)dW{t), (23) 

on (rj, 3", , Q). The initial condition is assumed to satisfy Assumption [3l 
Assumption □ is satisfied with d = l, K{t) = A- 'yX{t)B, R{t) = B. 

The key point of the construction of Section [2] and of its interpretation is 
the fact that is a martingale. To this end we compute its stochastic 

differential by using the Ito rules and we get 

dm)\m) = {<^m\m) + m^m) + {mmm) 

= {^{t)\ [A*+A- -fX{t) {B* + B) + B*B] ^{t))dt 

+ {m\{B* + B)m)dW{t). (24) 

Then, the process (||V'(OlP) can be a martingale only if the term in front of 
dt (the drift term) is equal to zero. This imposes that 

A* + A-jX{t){B* + B) + B*B = 0, Vt. (25) 

By taking the mean of this equation we get A* + A + B*B = 0; then, we 
need also B* + B. These conditions impose that there are two self-adjoint 
operators L and Hq such that B = —iL and A = —iHq — ^L"^. As a 
consequence the initial equation ([201) becomes 



dm 



-i{Ho-7X{t)L)-^L^ 



ij{t)dt - iLilj{t)dW{t) . (26) 



Now, being X{t) a continuous adapted process and Hq = Hq* G £(?C), 
L = L* e also Assumption [2] holds with 

H{t) = Ho--fX{t)L, R{t) = -iL, K{t) = -i{Ho--fX{t)L) -^L\ 



Moreover, we have Eq[|X(i)|] < ^EQ[X(t)^] < 1/^/27, 

Eq[\\H{t)\\] < \\Ho\\ + 7 ll^ll ^q[\Xm <\\Ho\\ + ^ \\L\\ , 
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which imphes condition (j3ap . Condition (jShj) and Assumption [J] are trivially 
satisfied because R{t) is non random, time independent and bounded. 

As all Assumptions [IHll hold, also all statements of Theorems [H [2] and 
Remarks [THl] hold. In particular the ISSE (p6]) has a pathwise unique solu- 
tion. 

What is peculiar of the present model is that Eqs. (p^ - ([26]) give 
= ||V'(0)||^ or that the probability densities are independent of time, p{t) = 
p{0), cf. Eqs. dS]) and ([H]), which give m(t) = and W{t) = W(t). We have 
also, from Eq. dSD, Tp{t) = ^(t)/ ||^(0)||, if \\ip{0)\\ / 0, and, from RemarkEl 
P\F) = q(F), for all events F G 3't, independent of Jq. As a consequence 
the change of probability has no effect (the new probability is equal to the 
initial for events independent of 3"o)- In other terms, no information has 
been extracted from the measurement interpretation. 

Moreover, the property ||V'(^)|| = llV'oll is in agreement with a purely 
Hamiltonian evolution. More precisely, let Texp{---} denotes the time 
ordered exponential; then, the formal solution of Eq. (I26p is given by 

^(t) = Texp|-i^ {Ho-jX{s)L)ds-iJ^ LdVF(s)|vo. 

The evolution of the quantum system is then completely determined by the 
time-dependent, random Hamiltonian 

Ht = Ho+ (W{t) - 7X(t)) L. 

Let us stress that it is a formal expression, due to the presence of W{t). 

This shows that the usual measurement interpretation of (j20p coloured 
with an Ornstein-Uhlenbeck process gives raise to a random Hamiltonian 
evolution. As announced, we recover the framework of the evolution of a 
closed system incorporating a random environment characterised in terms 
of an Ornstein-Uhlenbeck noise. 

One can investigate the evolution of the corresponding density matrices. 
To this end, we consider the pure state process {o'{t)) defined by (7{t) = 
\ijj{t)){'tp{t)\, Eq. d?]). By using Ito rules, the process satisfies the 

stochastic differential equation (SDE) 

da{t) = -i[Ho - -fX(t)L, a{t)]dt - i[L, a{t)]dW{t) - ^ [L, [L, a{t)]]dt, (27) 

which is, of course, equivalent to (j26p with random Liouville operator (jlOp 
given by 

C{t) = -i[Ho - ^X{t)L, ■]dt - \ [L, [L, ■]]dt. 
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Let us stress that the presence of the Ornstein-Uhlenbeck process imphes 
that the solution {cr{t)) of Eq. ([27|) is not a Markov process. 

Taking the expectation, we get the evolution of the mean r]{t) = EQ[a{t)], 
Eq. ([IS]), which turns out to be 

^r?(t) = -i[Ho,m - ^ [L,[L,7?(t)]] + i7[L,Eq[X(t)a(t)]]. (28) 

Note that it is not a closed master equation for the mean state r](t). Actually, 
we have derived a model with memory for the mean state. Indeed, the 
term i7[L, EQ[X(t)cr(t)]] introduces non-Markovian memory effects in the 
dynamics. Moreover, Eq. (|26p is an unravelling of the master equation (j28p . 

4 Projection techniques and closed master equa- 
tions with memory 

As we have seen in Eq. (|16|) . the a priori states or average states rj{t) = 
Eq[cr(t)] = Epf satisfy the equation ri{t) = Eq [£(t)[cr(t)]] , which is not 
closed because both C{t) and a{t) are random. However, at least heurist- 
ically, some kind of generalised master equations can be obtained by using 
the Nakajima-Zwanzig projection technique [21 Section 9.1.2]. 

Let us introduce the projection operators on the relevant part (the mean) 
and on the non relevant one: 

r[---]:=E^[---], Q:=t-V. 

Then, we have rj{t) = 'P[cr{t)] and we define the non relevant part of the 
state, the mean Liouville operator and the difference from the mean of the 
Liouville operator 

a^{t) := Q[ait)]=ait)-rj{t), 

Cuit) := Ec,[C{t)], AL{t) := L{t) - Cuit). 

By using the projection operators and the fact that the stochastic integrals 
have zero mean, which means V dWj{s) ■ ■ ■ = 0, from ([8]) we get the 
system of equations 

flit) = CuitMi)] + V o AC{t)[a^{t)], (29a) 
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da^it) = QoCit)[a^{t)]dt + Y^njit)[a^it)]dWjit) 

i=i 

d 

+ Q o C{tMt)]dt + J2'^j{tMt)]dWj{t). (29b) 
i=i 

As one can check by using Ito formula, the formal solution of Eq. ()29bp 
can be written as 

a^{t) = QoV{t,0)[a^m + QoV{t,s)o (^C{s) - Y^T^jisf^visMs 



+ QoV(t,0) 



fv{s,o)-'on,{sMs)]dw,{s) 



(30) 



where V{t, r) is the fundamental solution (or propagator) of the ISME ^ 
and satisfies the SDE 

V{t,r) = t+ dsC{s)oV{s,r) + Y dWj{s)TZj{s) o V{s,r). 

J r J r 

Let us stress that if one includes V(t, 0) into the stochastic integrals, from 
one side one gets the simpler expression V(t,0) o V(s,0)~^ = V{t,s). But 
the propagator is a stochastic process and could be non adapted. To over- 
come this difficulty one should use some definition of anticipating stochastic 
integral. So, we prefer the formulation with V{t, 0) outside the stochastic 
integral, in order to have only adapted integrands. 

By introducing the quantity pOj) into Eq. (|29ap . we get the generalised 
master equation for the a priori states 



rj{t) = Jit) + CuiiMt)] + / ICit,sMs)]ds 



where 



AC{t) o Qo V(t,0) 



Y / v{s,o)-'on,{sMs)]dWj{s) 



(31) 



J{t) := Eq [AC{t) o Q o V{t, 0)[a^{0)]] 

is an inhomogeneous term which disappears if the initial state is non random, 
i.e. (T_l(0) = 0, and 



JC{t,s) := Eq 



A£(t) oQoV{t,s) o (£(s) -^7^J(s 
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is an integral memory kernel. Also the last term in (j3ip is a memory con- 
tribution, linear in rj and depending on its whole trajectory up to t. 

Let us stress that to compute the terms appearing in Eq. (jSip and to 
solve it is not simpler than to solve Eq. ^ and to compute the mean of the 
solution. The meaning of Eq. (|3ip is theoretical: it is a quantum master 
equation with memory and Eq. ()14p gives an unravelling of it. 

While the best way to study a concrete model is to simulate the stochastic 
Schrodinger equation, Eq. ([3T|) could be the starting point for some ap- 
proximation. A possibility is to approximate V(t, r) by the deterministic 
evolution generated by the mean Liouville operator: 

VM(.t,r) = t+ [ dsCM{s)oVM{s,r). 

J r 

If we take also cr±(0) = 0, we get 



r7(t)~£M(i)[^(i)]+ / /Ci(t, s)[r/(s)]ds 

JO 

d ft 



+ E, 



A/:(t) 



V / VM(^,s)o7^,.(s)[r/(s)]dW^,■(s) 



(32) 



/Ci(t,s) := Ec 



A£(t)oVM(t,s)o (A£(s)-^A7^2(s; 



where A7^2(s) = TZjisf - EQ[nj{sf]. 

For the model of the previous section we have: TZjis) =1Z = —\[L, •], 

C-M{t) =Cm = -i[Ho, ■] - \[L, [L, •]], Vuit, s) = e^^^*"^). 



A£(t) 



■j=i 



^ / VM{t,s)on,{sMs)]dWj{s) 



X{t) [ dVF(s)7^oe^^(*-') o7^[r/(s)] 
Jo 

= -7/ ds7^oe(^M-T)(*-") o7^[^?(s)], 



/Ci(t, s) = 7=^ E(i[x{t)x{s)]n o e^M(i-s) o 7^ = 2 7^ o e(^^-T)(*-^) o 7^. 
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Finally, the approximation of the non Markovian master equation turns out 
to be 



+ ^ / dJL,e(^^-^)(*-^)[[L,r7(s)]]l . (33) 



However, we have no results on the positivity preserving properties of such 
an approximate evolution equation, while the unravelling of the complete 
equation guarantees complete positivity and feasibility of numerical simula- 
tions. 
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